Two Mode Photon Bunching Effect as Witness of Quantum Criticality in Circuit QED 
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We suggest a scheme to probe critical phenomena at a quantum phase transition (QPT) using 
the quantum correlation of two photonic modes simultaneously coupled to a critical system. As 
an experimentally accessible physical implementation, a circuit QED system is formed by a capac- 
itively coupled Josephson junction qubit array interacting with one superconducting transmission 
line resonator (TLR). It realizes an Ising chain in the transverse field (ICTF) which interacts with 
the two magnetic modes propagating in the TLR. We demonstrate that in the vicinity of criticality 
the originally independent fields tend to display photon bunching effects due to their interaction 
with the ICTF. Thus, the occurrence of the QPT is reflected by the quantum characteristics of the 
photonic fields. 

PACS numbers: 42.72.-g, 74.81.Fa, 73.43.Nq 



I. INTRODUCTION 



Inspired by the fast developments of quantum informa- 
tion [l|, 0] , quantum phase transition (QPT) Q has re- 
newed much attention in different fields of physics rang- 
ingfrom condensed matter physics to quantum optics 
[HIH . It was found that @ at the quantum critical point 
the dynamic evolution of a quantum critical system is 
so extremely sensitive that it can enhance the quantum 
decoherence of an external system coupled to it. This 
ultra-sensitivity is characterized by the Loschmidt echo, 
which is a well-known concept in quantum chaos [7[. In 
this sense, the quantum-classical transition from a pure 
state to a mixed one is induced by the quantum criticality 
of this surrounding system. This discovery motivated a 
new scheme to probe the QPT by exploring the quantum 
coherence in the external system and its losses fl2l ]. 

Moreover, this probing mechanism for quantum crit- 
icality was illustrated by a circuit QED architecture 
d, M, which was formed by a superconducting 

Josephson junction qubit array interacting with a one- 
dimensional superconducting transmission line resonator 
(TLR) [12]. The superconducting qubit array was mod- 
eled as an Ising chain in the transverse field (ICTF). This 
investigation showed that the QPT phenomenon in the 
superconducting qubit array was evidently revealed by 
the correlation spectrum of TLR output. Though this 
mechanism for the circuit QED system has not been ex- 
perimentally tested, an NMR simulation experiment [l3[ 
has been carried out to demonstrate the QPT-likephe- 
nomenon (energy level crossing ) as predicted in ref. [6| by 
exploring the increased sensibility of the quantum system 
to perturbation when it is close to a critical point. 

For the above circuit QED architecture to demonstrate 
the probing of the QPT, we notice that with two modes 
simultaneously coupled to a charge qubit, their squeez- 
ing effect was investigated theoretically [ll|. Here, we 
consider the full application of quantum optics approach 
[rij in the detection of the QPT by considering the higher 



order quantum coherence. To this end, we consider that 
a Josephson junction qubit array modeled as the ICTF 
simultaneously couples to two modes propagating in the 
TLR. Since all quasi-spins homogeneously interact with 
the fields, we can obtain the first (second) order correla- 
tion function of the two fields. According to our calcu- 
lation, the second order quantum coherence is given in 
terms of the decoherence factor of the ICTF. As proven 
in the Appendix, the norm of the decoherence factor de- 
creases exceptionally when the ICTF is at the critical 
point. Therefore, the photon bunching effect occurs since 
the second order quantum coherence of the steady state 
is smaller compared with its initial value. And these re- 
sults show genetic characteristics of the quantum spin 
chain in the vicinity of the critical point. 

The paper is structured as follows. The next section 
describes the ICTF formed by a capacitively coupled 
Josephson junction qubit array is coupled to two inde- 
pendent fields propagating in the TLR. Then the detec- 
tion scheme of the QPT and the correlation functions of 
two mode fields are given in Sec. III. A brief summary 
is concluded in Sec. IV. Furthermore, in addition to the 
main body of the paper, Appendix [A] presents details 
about the calculation of the decoherence factor. 



II. CIRCUIT QED BASED SETUP FOR 
PROBING QUANTUN CRITICALITY 

We consider a circuit QED system illustrated in Fig.l. 
N Cooper pair boxes (CPBs) are capacitively coupled 
one by one. Formed by a superconducting island con- 
nected with two Josephson junctions, each CPB is a di- 
rect current superconducting quantum interference de- 
vice (dcSQUID). Since the magnetic flux $> x threading 
the dcSQUID is tunable, the effective Josephson tun- 
nelling energy can be varied. With proper bias voltage, 
the CPB behaves as a qubit near the degeneracy point 
and then Josephson junction qubit array becomes a spin 
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FIG. 1: (color online) Schematic diagram of superconducting 
Ising chain interacting with two largely detuned modes (u>i = 
3u>2) individually transmitted in the TLR. Since the CPBs 
are located at the antinodes of both modes, i.e, Xj = (j — 
1/2)ttc/uj2, they only interact with the magnetic fields. 



chain with N 1/2-spms. When the coupling capacitance 
C m between two CPBs is much smaller than the total 
one Cs to each CPB, the high order terms in Hamil- 
tonian can be neglected and only the nearest neighbor 
interaction is considered. Then the qubit array can be 
approximated as an ICTF with [l2| 



N 



'3 U 3 + 



l). 



(1) 



where o x = -|0)(1| - |1)(0| and a* = |0)(0| - |1)(1| 
with \n) being the state of n extra Cooper pair on the 
superconducting island, A = B x /B and B = e 2 C m /C^, 
B x — Ej cos(ir$ x /&o) is the Josephson energy of each 
CPB with Ej the Josephson energy of single junction, 
$0 = h/2e the flux quantum. 

In a one dimensional TLR, the electric current and 
voltage at the position x are given as 



/(M) = £V^(4 + ^sm(^), (2) 



V(x,t) 



\l —{al-a k )cos{— ), 



(3) 



where at is the creation operator with frequency uik, L 
the length of TLR, I and c the induction and capacitance 
of per unit length of TLR respectively, k positive integer. 
Therefore, a CPB located at the antinode is only coupled 
to the magnetic field since the electric field vanishes. Ac- 
cording to Ampere's circuital law, when a dc SQUID loop 
is placed at a distance r with respect to the center of the 
TLR, the quantum magnetic flux that threads it is 
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Ojfe)sin(-g-), 



(4) 



where fio is the vacuum magnetic permeability, 5* the 
area of dc SQUID loop. 



The interaction between the CPBs and the magnetic 
field is written as 



C0S(— 1^)<7 



(5) 



In our consideration, two independent modes with fre- 
quencies u)\ = 3co>2 are propagating in the TLR. All the 
CPBs are placed the antinodes of the both modes with 
the positions 



(j - 1/2)™ 



t0 2 



where 



U)2 



TTVL 

IvP 



(6) 



(7) 



j and M are positive integers, v the velocity of the light. 

Since <I> g <C $0i under the rotating wave approxima- 
tion [l4| , the interaction Hamiltonian is approximated to 
the second order, 

H^=Ej[{l-^ n l-^ n l)-r } la\a 1 -r ) lala 2 ]Y J ^h ( 8 ) 



where the coupling constants between the two modes and 
individual spins are 



TT/IqS / Kul 



27rr$ V LI 



(9) 



For realistic parameters, Ce = 600 aF, C m — 30 aF, 
L Q = 1 cm, S = 10^m 2 , r = 1/xm, N = 500, Ej = 13 
GHz, we have B = 1.6 GHz, u 2 « 120 GHz, 771 = V3?72, 
and J] 2 « 0.1 [3. 

Thus, the total Hamiltonian is written as 

H =Ej[(l - )-rj{ - ]-<nl) - r)la\ai - "<]\a\a 2 \Y^o x 3 



e 2 G, 



r 2 



" VjVj+i + wiojoi + uJ2a\a-2- (10) 



Furthermore, since there is no energy exchange be- 
tween the fields and the ICTF, the total Hamiltonian 
can be decomposed into invariant subspaces with respect 
to the Fock state of the fields, 



H = ^H {m ' n) \m)\n){n\(m\ 



(11) 



where 



N 



H {m,n) = B J2(\ m , n <7* + a|<7? +1 ) (12) 

i=i 

with B = e 2 C m /Cl and A m „ = E ,[l - (m + 1/2)77? - 
(n + 1/2) V 2 )]/B. 
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Generally speaking, the Hamiltonian of ICTF Ho is 
transformed into a quadratic fermion form with Jordan- 
Wigner transformation 0] 



fc=i 

Then, by introducing quasi-particle operator [l| 

N 



Ik 



E 



[cos y Cj ■ - z sin — cj] 



AT 



with 



#fc(A) = tan 



sin fc 



A — cos k ' 
Ho is diagonalized as 

fc 

with single particle spectrum being 



e fc (A) = 2Bv / l + A 2 -2Acosfc. 



(13) 



(14) 



(15) 



(16) 



(17) 



And the ground state \G) corresponds to no quasi- 
particle excitation at all. 



III. PHOTON BUNCHING EFFECT 

Followed by a series of advances, i.e., resonance fluores- 
cence, the Hanbury-Brown-Twiss experiment fltl ] reopen 
philosophical debate about photons [I?) and set itself as 



the milestone in the development of quantum optics. All 
these experimental phenomena are associated with the 
correlation functions of the field. Here, we consider it 
as the method to detect the QPT since the two fields 
propagating in the TLR interact with the quasi-spins re- 
spectively. 

First of all, we define an operator 



A = a± + ia 2 . 



(18) 



The first order correlation function is written as 
(A* (t)A). Here, the bracket (■■•) denotes average over 
the initial state, with the Ising chain in the ground state 
\G) and the two fields being in arbitrary pure states 
^ m c m |m) and ^2 n d n \n) respectively. Therefore, 



(A\t)A)=J2\c m \ 2 \d n f(mr^ 



(m,n) 



i) 



- ict n _ 1 c m (% l+1 d n y/m(n + l)r^" ^; n+1) 

+ * c m+l c ™eC-l<W ( m + ^ r m,n-T 

(19) 



where 



r%':)(t) = {G\e iH ^ n %- 



iff* 



h \G) (20) 



is the decoherence factor Q which measures the overlap 
of the ground state evolving under two different Hamil- 
tonians. Details about its calculation is presented in Ap- 
pendix[A] In Ref . [l2| . it was discovered that for the same 
amount of environment dissipation the first order corre- 
lation function of the single mode decreased more rapidly 
in the vicinity of the QPT than in the other region. 

Moreover, the second order correlation function 
(A' A* (t)A(t)A) is analytically written as 



E 

m.n 



{rf A\t)A{t)A) = 

\c m \ 2 \d n \ 2 [(m + n)(m + n - 1) + mnr^f e~^-^ + mnr^e^-" 1 '*] 

+ ic^Cmd^Mim + n- l)v/(m + l)n + [my/(m + l)^™^ + (n - 1) y/{m + l)nr^f- V^ 1 "^*} 

- icl_ lCm d* n+1 d n {(m + n- l)y/m(n + l) + [(m 1) y/m(n + l)^" + n y/m(n + l>^f ^e^ 1 "" ^} 
~ c* m+2 c m d* n _ 2 d n ^{m + 2)(m + l)n(n - 1)r ^+^) e ^^)t 

- c* m „ 2 c m d* n+2 dW(™ l)m(n + l)(n + 'e^^)' (21) 



r 



Thus, for the fields initially in the state (|0) + |l))/\/2, 
it is straightforward to obtain 

(A^ A\t)A{t)A) = ~[1 + Re(r£i ( V^ 1 -" ;2 W)], (22) 

where Re(x) means the real part of x. 



As proven in Appendix \X[ in the vicinity of the QPT, 
the square of the norm of r^f* 1 (t) decreases more rapidly 
than exponential, i.e, 



|rfe 0) (*)| 2 <e 



(23) 



where 7 = 4B 2 (A li0 - A 04 ) 2 £;(fc c )/(Ao,i - l) 2 , E(k c ) 
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FIG. 2: (color online) The decoherence factor \r^f (t)\ 2 for 
both fields in (|0) + \1))/V2 is plotted with N = 8000. Blue 
dashed line for A = 1, red solid line for A = 0.1, and green 
dotted line for A = 2. In all figures, t is in units of 1/B. 



Att 2 N c (N c + 1)(2N C + l)/6N 2 with N c being the nearest 
integer to Nk c /2-K. 

It can be seen that there is a vanishing numerator 
E(k c ) as N — > oo. It is doubtful that the exponential 



decay of 



.CM) 



(t)\ 2 can truly occur since the QPT takes 



place in the thcrmodynamical limit. However, as the size 
of the ICTF gets larger, we can adjust the parameter 
Ao.i closer to the critical point to make the denomina- 
tor (Arj.i — l) 2 small enough. In that case, 7 stays as a 

constant and the Irg 1 ^ ' (t)\ 2 decreases exponentially with 
time. For a real system, N is finite for the demonstration 
of the QPT. To test the validity of the above analysis, 
we resort to numerical simulation. In FigJ^l we plot the 

evolution of |f'Q 1 j ' ) (£)| 2 according to Eq. (|A9p . It can be 

„(i.o) 



seen that despite some oscillations |tq ± i'(t)\ 2 decays ex- 
ceptionally at the critical point. 

According to Ref. [3, the photon bunching and anti- 
bunching effects are associated with the second order de- 
gree of coherence 



9 {2 \t)- 



{A*At{t)A(t)A) 
{A\A){A^t)A{t)Y 



(24) 



which is the normalized second order correlation function 
of the fields. For the fields both in the state ( 1 0) + 1 1) )/ y/2, 
the second order degree of coherence is simplified as 



(2) W = - 



1 1 + Re(r { i 



(1.0) i{u!-u%)t 



1 



Im(r^{ 0)ei( " 1_W2) *) 



(25) 



with Im(x) being the image part of x. 

Since the norm of the decoherence factor decreases ex- 
ponentially at the critical point, it is obvious that both 
the real and image parts of r ^-^> e ^i-^)t will vanish in 
that limit. As a consequence, we expect the second order 
degree of coherence to be less than unity in the steady 
state, i.e., g {2 \t) = 1/2 < # (2) (U) = 1. Generally speak- 
ing, classical fields such as thermal light and coherent 




10 20 30 40 50 60 70 80 90 100 
1 

(b) 

I I 1 1 1 1 1 1 1 1 1 1 



10 20 30 40 50 60 70 80 90 100 
t 

(c) 



0.01 




10 20 30 40 50 60 70 



FIG. 3: (color online) The second order degree of coherence 
g (2) (t) for N = 4000 and (|0) + \1))/V2 is plotted with (a) 
A = 1, (b) A = 0.1, (c) A = 2. 



light, prefer to distribute themselves in bunches rather 
than at random. They exhibit less correlation for time 
longer than the correlation time. This is so called bunch- 
ing effect [3]. On the contrary, in certain quantum op- 
tical systems, fewer quantum photons are detected close 
together than further apart. And the photon antibunch- 
ing observed in fluorescent light from a two-level atom 
171 ] is of such kind. Here, since the two fields involved 
are two independent modes, we expect photons to be nei- 
ther bunching nor antibunching, regardless of quantum 
mechanical fields or classical fields. However, as shown 
in FigJ3J when the Ising chain is at the critical point, 
the two independent fields initially in (|0) + |l))/\/2 dis- 
play the photon bunching effect. Further witness is also 
demonstrated in FigfJJa-c). It can also be proven that 
g^ 2 \t) < g( 2 \Q) for both fields in the coherent state \a) 
which is not shown here. In FigfJJd), we plot the time 
evolution of the second order of coherence for this case. 
Here, we remark that the two initially independent quan- 
tum fields display the classical effect due to their com- 
mon interaction with the quantum critical system. As 
illustrated in Eq. (|20|) , two initially identical states evolve 
under two slightly different Hamiltonians. Although the 
difference between these Hamiltonians are tiny, their evo- 
lution trajectories are quite distinct in the vicinity of the 
QPT. Thus, this slight difference leads to the exponential 
decay of their decoherence factor. It can be understood 
as a signature of quantum chaos jf|. 

Furthermore, for the parameters mentioned after 
Eq.© and N c — N/10, both the real and imaginary 
parts of r ^ 1 ^°) e , ( LJ i- w 2)t decay with a rate of the order 
*/7 ~ 2.5 GHz. Since the dissipation rate of the first 
excitation mode is about 6.3 MHz Q, we can neglect the 
influence due to the dissipation of TLR. 
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FIG. 4: (color online) The second order degree of coherence 
g {2) (t) is plotted at the critical point for (|0) + |1»/V2 with 
(a) N = 2000, (b) N = 4000, (c) N = 8000. For (d), both 
fields are in the coherent states ja) with a = 1 and N = 8000. 
Note that at the steady state g^ 2 \t) is a little smaller than 
its original value 1 as indicated by the red horizontal line. 



IV. CONCLUSION AND REMARK 

To conclude, we have explored the possibility to probe 
quantum criticality in the ICTF by detecting the higher 
order quantum coherence of the two modes of cavity fields 
coupled to the spins. We suggest a physical implemen- 
tation of this theoretical scheme based on a circuit QED 
system where the capacitively coupled CPBs are coupled 
to the TLR. Situated at the antinodes of both modes 
propagating in the TLR, CPBs are only coupled to the 
magnetic fields. In a heuristic way, we show the decoher- 
ence factor decays exponentially with time in the vicinity 
of the critical point. The second order of coherence is 
smaller than one at the steady state. Thus, the two ini- 
tially independent modes demonstrate photon bunching 
effect. This can serve as a witness of the QPT. 

On the other hand, we have not investigated decoher- 
ence originated from the dissipation of the CPBs. We 
notice that in a recent work [18(], the QPT in the dis- 
sipative random transverse-field Ising chain was investi- 
gated. It was discovered that the quantum critical point 
was ruined by the interplay between quantum fluctua- 
tions and Ohmic dissipation. Further exploration may 
be done when such kind of effect is considered. 
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APPENDIX A: DECOHERENCE FACTOR 

Following the method introduced in Ref. [l2j , the deco- 

herence factor r^?'"? can be calculated in the following 
way. 

By introducing the spin-1 pseudospin operators [l9j 

Sxk = i{l-klk+ll k lt), 

Syk=ll k ll-l-klk, 

Szk = lllk + l- k l-k ~ 1, (Al) 

the Hamiltonian Ho can also be rewritten as 

H = ^e k s zk . (A2) 

fc>0 

Because there is no energy exchange between the 
two modes and the qubit array, the total Hamilto- 
nian can be decomposed into invariant subspaces with 
respect to the Fock state of the fields, i.e., H = 
Y, m , n H {m ' n) \m)\n){n\{rn\, where 

N 

H (m,n) = B J2(\ m ,n<7] + (A3) 
3=1 

With the pseudospin operators, we can also diagonalize 
the Hamiltonian as 



(m,n) (m,n) 
k > 



(A4) 



k>0 



2a 



(m,n) 
zk 

(m,n) /j(m,n) 



where s,u''"' = s zk cos2al 



(m,n) . n (m,n) ■ . , 

r k + s xk sin 2a k , with 

g Jm,n) _ _ ( \ \ /)(m,n) _ 
'ki £ k — £k\^m,n), v k — 

0k(^m,n)- 

Therefore, the ground state of Hq is the product state 
of all pseudospins down \—}k, 



i G >=n 

fc>0 

= IK 



cosa l ' |->* 



sin a 



(m.n) I \ (m.n) 





k>0 










with |±) 


(m.n) 


being 


the ei; 


'en states of 


(m,n) 
s z)c 


Since 








(m,n) 
k 


(±l±) 


(m',n') 
k> 


= S k k' 


/ (m.n) 

cos(a^, ' - 




(m,n) 
k 




(m'.n') 

y 


= Skk' 


/ (m,n) 

sm(a k — 




(m,n) 
k 




(m',n') 
k' 


= Skk' 


. / (m ,n ) 

sm(<4 


-4 m < n) ), 



) 

(A5) 
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For heuristic analysis, wc obtain the short time behav- 
ior of |ToY^(*)| 2 at tne critical point. 



with 


















°+,+,fe 


,n') 


(m. 

= sin ajj. 


n) 


(m 

cosaj. 


,n') ■ i 

sm( 


(m,n) 






°-,-,fc 


,n') 


= cos ajj. 


,n) 


fm' 

sma k 


,n') ■ i 

sm( 


a <V,n') _ 


-4 m,n) ; 


°+ 


,n') 


(m. 

= sin ' 


n) 


(m f 
sma fc 


cos( 


a ( m ,„) _ 


(m' 




°-,+,fc 


,n') 


= cos a k 


,») 


cos a]. 


>') 

cos 


/ (m,n) 


-a* 





we have 

«w = n E cixr' B ' ) e <(o *' im,B,+ *-' eim ''"' ,)t 

fe>0 afe,6fc=± 

(A7) 



k>0 

= n [^ a (4 1,0) - 4 0,1) ) cos(4 M) + 4 0,1) )* + 4 0,1) ) cos( £ ^ 0) - 4 0,1) )*] 2 

fe>0 

+ [sin( a ^°) + 4 ' 1 ) ) *n(4 1,0) " 4 M) ) cos^ + e™)t 

- cosfa™ + 4 0>1) ) costal - a?' 1 ') an^ + ef 1 ^] 2 (A9) 



Since all factors Fk of \tq'^ (t)\ 2 have a norm less than Since 
unity, we may expect the |7n\' ^(t)| 2 to vanish under cer- 

' (1 0) (0 1) (1 0) (0 1) (1 0) 

tain conditions. Here, we set a cutoff frequency k c and w e fe ' — e fc ' «£ t ! + e fc ' w 2e fc ' , 
hence wc have \4\°\t)\ 2 < fl For sma11 fc , 

we have 

fc>o wc focus on the short time behavior and therefore 



.(i,o) 
fc fc 

_(o,i) 
£ k 

e k (\) 



,,(0.!) 



2B|1-A ll0 |, Ir^WI 2 
2B|1- Ao,i|, 



fee 



< n [i - 2 (4 1,0) - 4 0,1) ) 2 ] + 2(4 M) - 4 0,1) ) 2 cos(24 1 '°)i ) 



A-l' fc>o 
V X^-l - A^< = ft 1 2(4 1,0) - 4 ' 1} ) 2 S in 2 (4 1,0) *) 



I, k k 



k>0 



2^Ao,i-l Ao,o-l^ -rv, n fc 2 (Ai,o - A 1) 2 2/ OR+ i\ m 

= ll 1 ~ 2 ( Alo _ 1 ) 2 ( Aoi _ 1 )2 sm (2flf|Ai, -l|). 
To the second order of ai}' ^ ± a[°' 1 ' ) , we obtain fe>0 



l r o,i (*)l As Ai, -> 1, we have 

< ft [1 - 2(«r - 4°' 1) ) 2 ] cos 2 (4 1 '° ) - 4 M) )* \ r ™ (t) , a < e-* 2 , (A10) 



fe>0 



+ 2K -« fc ) cos( £fc + £fc )icos(e fc -e fc )i where 7 = AB 2 {Xl Q _ Xq ^2^/^ i _ 1)2 ^ E{kc) = 
+ [1 - (4 1 ' )) 2 - (4° 4) ) 2 ] sin 2 (4 ll0) - 4° 4) ) i 4tt 2 A^ c + 1)(2A C + 1)/6A 2 with A^ c being the nearest 



-2[(4 1 -°)) a -(ar) a ]™(4 1,0) -4 0,1) )' 

xsin(4 1 '°) + 4°' 1 )) t . 



integer to Nk c /2n. 
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